MovingPandas Minimum Viable Example

MovingPandas provides a trajectory datatype based on GeoPandas. The project home is at https://github.com/movingpandas/movingpandas

The documentation is available at https://movingpandas.readthedocs.io/en/main/

import warnings

warnings.filterwarnings("ignore")

import pandas as pd
import geopandas as gpd
import movingpandas as mpd
from hvplot import pandas
from datetime import datetime, timedelta
import folium

Loading trajectory data from a GeoPackage

The MovingPandas repository contains a demo GeoPackage file that can be loaded as follows:

gdf = gpd.read_file("../data/geolife_small.gpkg")
gdf.plot(figsize=(9, 5))

<Axes: >

gdf.hvplot(geo=True).opts(active_tools=["pan", "wheel_zoom"])

c:\Users\Agarkovam\AppData\Local\miniforge3\envs\mpd-ex\lib\site-packages\geoviews\operation\__init__.py:14: HoloviewsDeprecationWarning: 'ResamplingOperation' is deprecated and will be removed in version 1.18, use 'ResampleOperation2D' instead.
  from holoviews.operation.datashader import (

gdf.explore()

Make this Notebook Trusted to load map: File -> Trust Notebook

After reading the trajectory point data from file, we want to construct the trajectories.

Creating a TrajectoryCollection

tc = mpd.TrajectoryCollection(gdf, "trajectory_id", t="t")
tc

TrajectoryCollection with 5 trajectories

tc.plot(column="trajectory_id", legend=True, figsize=(9, 5))

<Axes: >

tc.explore(column="trajectory_id", cmap="plasma", style_kwds={"weight": 4})

Make this Notebook Trusted to load map: File -> Trust Notebook

Exploring movement speed

tc.plot(
    column="speed", linewidth=5, capstyle="round", legend=True, vmax=20, figsize=(9, 5)
)

<Axes: >

Detecting stops

detector = mpd.TrajectoryStopDetector(tc)
stop_points = detector.get_stop_points(
    min_duration=timedelta(seconds=120), max_diameter=100
)

ax = tc.plot(figsize=(9, 5))
stop_points.plot(ax=ax, color="red", markersize=100)

<Axes: >

(
    tc.hvplot(line_width=7, tiles=None, frame_width=400, frame_height=400)
    * stop_points.hvplot(geo=True, color="black", size=100)
)

m = tc.explore(
    column="trajectory_id",
    cmap="autumn",
    style_kwds={"weight": 4},
    name="Trajectories",
)

stop_points.explore(
    m=m,
    marker_kwds={"radius": 4},
    name="Stop points",
)

folium.TileLayer("CartoDB dark_matter").add_to(m)
folium.LayerControl().add_to(m)

m

Make this Notebook Trusted to load map: File -> Trust Notebook